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ATOMIC  APPROACHES  TO  DEFECT  THERMOCHEMISTRY 


J.A.  Van  Vechten  and  J.F.  Wager 
Oregon  State  University 

GOAL 

To  apply  atomistic  thermodynamic  theory,  Monte  Carlo  simulation,  and  experimental  analysis  to 
elucidate  the  identity  of  point  defects  in  semiconductors  as  well  as  to  understand  their  static  and  dynamic 
properties. 


TECHNICAL  ACCOMPLISHMENTS 

1.  Direct  Simulation  of  Atomic  Diffusion  Including  the  Effects  of  Carrier  Injection  and  Junction  Mis-Fit 
Strain:  "Resolving  Defect  Mega-Controversy  by  Glga  Event  Monte  Carlo  Simulation  of  the 
Macro-Consequences  of  Atom-Level  Assumptions  on  Microcomputers  at  Nano-Cost,"  J.A.  Van 
Vechten,  U.  Schmid,  and  N.C.  Myers  in  Defect  Control  in  Semiconductors,  ed.  by  K.  Sumino  (North 
Holland.  Amsterdam  1990)  p.  41;  “VIDSIM  -  A  Monte  Carlo  Program  for  the  Simulation  of  Atomic 
Diffusion  In  Diamond  and  Zinc-Blende  Structures,"  U.  Schmid,  N.C.  Myers,  and  J.A.  Van  Vechten, 
Comp.  Phys.  Comm.  58,  329  (1990);  "Failure  of  the  “Kick-Out"  Model  for  the  Diffusion  of  Au  into  Si 
When  Tested  by  Monte  Carlo  Simulation,"  U.  Schmid,  J.A.  Van  Vechten,  N.C.  Myers,  and  U.  Koch, 
Mat.  Res.  Symp.  Proc.  163,  609  (1990);  "Atomic  Diffusion  with  Strain  and  Injection,"  J.A.  Van 
Vechten,  Mat.  Res.  Soc.  Symp.  Proc.  184,  165  (1990);  "Surface  Treatment  Effects  on  Atomic 
Diffusion  in  Si  Explained  Without  Self  Interstitials,"  J.A.  Van  Vechten,  U.  Schmid,  and  Zhang  Q.-S., 
J.  Electron.  Mater  20,  431  (1991). 

A  major  and  long  term  goal  of  this  effort  is  to  establish  a  rigorous  understanding  of  the  relationship 
between  atom  scale  mate'ial  properties  (such  as  the  strengths  of  bonds  between  various  types  of  atoms, 
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vacancy  and  interstitial  formation  and  migrations  energies  energies)  and  device  scale  properties  that 
determine  the  performance,  reliability,  and  yield  of  electronic  devices  as  well  as  the  parameters  required  for 
their  production.  There  is  a  richly  developed  body  of  experimental  techniques  and  a  plethora  of  theoretical 
approaches  to  deal  with  the  atom  scale  class  of  properties  in  semiconductors  and  other  electronic  materials, 
but  only  a  very  tenuous  bridge  from  them  to  the  device  scale  problems.  Device  scale  problems  are 
inherently  complex.  Because  of  this  complexity  and  of  the  frantic  pressures  of  economic  competition,  device 
development  has  tended  to  fail  to  make  use  of  the  atom  scale  information  and  instead  to  rely  on  cut-and-try 
approaches. 

The  cut-and-try  atmosphere  of  device  development  has  spawned  a  number  of  semi-empirical  models 
for  the  critical  diffusion  and  interaction  processes  which  are  parameterized  to  fit  the  data  base  at  hand. 
They  assert  formulae  with  no  rigorous  derivation  and  assign  atom  level  parameters  with  no  regard  to  the 
results  of  atom  level  experiments.  Most  of  these  fail  badly  when  extrapolated  to  the  lower  temperatures, 
smaller  dimensions,  non-one-dimensional  diffusions,  and  much  sharper  gradients  characteristic  of  advancing 
technology.  They  also  tend  to  contradict  one  another  and  promote  unproductive  bickering. 

We  noticed  that  the  advance  of  computing  technology,  particularly  of  microcomputer  power  and  availability, 
make  it  now  practical  to  do  direct  simulations  of  atom  scale  processes  on  the  scale  of  modern  devices.  As 
several  billion  calculations  are  required  for  each  interesting  case,  such  a  project  was  heretofore  regarded 
as  impossibly  tedious.  We  have  developed,  demonstrated,  and  published  the  VIDSIM  program  which 
bridges  the  gap  between  the  atom  level  and  the  device  level  on  microcomputers  at  very  low  incremental 
cost  and  in  reasonable  time.  Unlike  almost  all  other  Monte  Carlo  simulations,  VIDSIM  does  not  invoke  any 
differential  equations.  It  simply  works  out  the  macroscopic  consequences  of  whatever  set  of  atom  level 
assumptions  and  parameters  the  user  wishes  to  use.  We  submit  that  most  of  the  bickering,  failure  of 
extrapolations,  and  wasted  effort  in  the  field  should  be  eliminated  when  all  concerned  can  see  for  themselves 
what  does  and  what  does  not  rigorously  follow  *rom  the  atom  level  io  thw  device  levs*  or  vice  versa. 

Our  first  major  success  in  this  effort  has  been  the  elimination  of  all  credibility  from  the  so  called 
"kick-out“  mechanism  for  atomic  diffusion  in  Si.  This  model  assumes  that  host  interstitials  mediate  the  major 
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part  of  host  and  dopant  diffusion  by  exchanging  places  with  host  or  impurity  atoms  on  substitutional  sites. 
Advocates  usually  assume  that  vacancies  are  not  able  to  annihilate  these  host  interstitials.  We  have  shown 
that,  contrary  to  wide  spread  belief  based  on  faulty  analysis  of  differential  equations,  the  “kick  out“ 
mechanism  of  diffusion,  can  not  produce  a  two  sided,  "U  shaped",  profile  from  a  one  sided  in-diffusion  of 
Au  or  other  transition  metals  that  diffuse  primarily  as  impurity  interstitials  but  occupy  predominately 
substitutional  sites.  We  have  identified  the  error  in  the  previous  analysis  (of  others)  which  led  to  the  false 
conclusion.  We  have  devised  a  computer  “bug",  which  when  introduced  to  the  simulation  can  force  the 
inconsistent  assumptions  of  the  previous  analysis  to  occur  (in  a  plair'y  unphysical  manner)  which  will  then 
produce  the  result  claimed  by  these  proponents  where  the  unphysical  condition  holds.  As  two  sided  profiles 
are  firmly  established  in  certain  cases,  this  result  eliminates  a  major  class  of  models  that  have  been 
advocated  to  explain  them. 

We  have  shown  that  the  fact  that  one  sided  in-d!ffuslon  of  Au  does  produce  a  two  sided  profile  of 
substitutional  Au  does  prove  that  SI  self  interstitials  can  play  no  significant  role  in  thermal  processes  in  Si. 
Furthermore,  several  surface  treatment  effects  which  others  have  claimed  require  Si  self  Interstitials  can  be 
naturally  explained  with  out  them.  (See  J.  Electron.  Mat.  20,  431  (1991),  which  Is  in  a  special  issue  devoted 
to  diffusion  and  which  was  edited  by  A.F.W.  Willoughby,  an  acknowledged  expert  on  the  subject.) 

The  kick-out  mechanism  has  also  been  invoked  for  diffusion  in  GaAs  and  other  lll-V  semiconductors. 
We  are  working  to  resolve  that  issue  now. 

We  have  shown  that  a  simple  version  of  the  Frank  Turnbull  mechanism  for  atomic  diffusion  does 
produce  a  two  side  profile  of  substitutional  impurities  from  a  one  sided  source,  as  is  often  observed. 
Moreover,  the  detailed  shape  of  the  profile  depends  on  surface  boundary  conditions  in  a  degree  and  manner 
not  previously  appreciated  in  the  field.  By  adjusting  assumptions  In  a  plausible  way,  the  simple  model  can 
be  made  to  fit  experimental  data. 

We  have  made  what  seem  to  be  first  ever  observations  of  the  distinctions  between  the  result  of 
In-diffusion  from  an  absorbing  boundary  condition,  which  seems  to  be  the  physically  realist,c  case  for  on 
Interface  to  a  second  phase  (e.g.,  vacuum  or  a  melt  or  a  eutectic),  and  the  coherent  boundary  to  a  surface 
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of  constant  concentration.  The  latter  was  the  case  considered  by  Einstein  in  his  study  of  Brownian  motion. 
Weil  known  results  of  Einstein's  study,  e.g.,  the  error  function  complement  relation,  do  not  occur  for  the 
absorbing  boundary  condition,  but  the  distinctions  are  subtle  if  one  does  not  know  a  priori  the  magnitude 
of  the  diffusivity  of  the  diffusing  species. 

2.  Independent  Measurement  of  Vacancy  Concentration  and  Diffusfvrtv  bv  Precise  Time-Dependent 

Measurement  of  Thermal  Expansion  at  Constant  Temperature:  R.K.  Graupner,  T.K.  Monson  and  J.A. 

Van  Vechten,  TBP. 

Advances  In  interferometry  over  the  past  decade  have  made  it  practical  to  measure  variations  in  the 
dimension  of  a  macroscopic  sample  with  a  precision  better  than  one  Angstrom  mean  surface  displacement 
over  periods  of  many  hours.  Thus,  one  can  now  measure  the  expansion  or  contraction  of  a  single  crystal 
sample  held  at  constant  temperature  as  vacancies  (and  interstitials  if  significant)  diffuse  in  or  out  of  the  bulk 
to  equilibrate  the  concentrations  of  these  point  defects  with  their  source/sink,  which  is  the  sample  surface. 
By  measuring  the  absolute  value  of  this  expansion  at  a  particular  time,  one  can  deduce  the  total  number 
of  vacancies  (minus  interstitials)  present  in  the  sample  at  that  time.  By  measuring  this  as  a  function  of  time, 
one  can  extract  the  diffusivity  of  these  native  point  defects  independent  of  their  concentrations.  (This  is 
regarded  as  a  major  coup  in  the  field  because  previous  measurements  give  only  the  product  of 
concentration  times  diffusivity  and  some  ambiguity  as  to  whether  vacancies  or  self  interstitials  are  dominant.) 
By  doing  these  measurements  at  different  temperatures,  one  can  detect  the  presence  of  multiple 
components  (e.g.,  interstitials  as  well  as  vacancies). 

Because  the  concentration  of  vacancies  expected  in  Si  and  other  tetrahedral  semiconductors  is  in 
the  10  parts  per  million  range  or  lower  at  manageable  temperatures,  there  would  at  first  seem  to  be  a  major 
problem  maintaining  T  so  constant  that  thermal  variation  of  the  lattice  constant  as  T  drifts  does  not  mask 
the  effect  one  is  seeking.  A  crude  estimate  implies  that  for  1  cm  sized  samples,  one  would  have  to  hold 
T  constant  to  0.01  C  at  1000  C  for  20  hours.  This  would  be  quite  challenging.  We  have  devised  a  way 
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around  this  problem  which  offers  also  an  opportunity  to  advance  the  art  of  non-invasive  semiconductor 
temperature  metrology.  We  cut  multiple  samples  each  In  their  equilibrium  shape  and  the  same  size  from 
one  single  crystal  ingot  of  the  material.  We  have  devised  a  two  zone  furnace  with  vacuum  transfer 
mechanism.  We  will  place  one  of  the  sample  in  each  of  zones  and  leave  them  there  long  enough  to 
eauilibrate  at  their  respective  temperatures.  We  then  transfer  the  sample  under  test  repeatedly  between  the 
two  zones.  We  measure  the  variation  in  the  dimension  of  the  sample  under  test  with  respect  to  the 
equilibrated  sample  in  Its  zone.  In  this  way  any  drift  In  T,  or  net  evaporation,  etc.,  affects  both  samples 
similarly  and  we  measure  only  the  effect  of  the  equilibration  of  the  sample  under  test. 

In  doing  this,  we  will  have  the  option  of  applying  feed-back  to  the  furnace  heating  system  so  that 
the  dimension  of  the  equilibrated  sample  is  held  constant  to  very  high  accuracy.  This  will  allow  us  to 
prepare  samples  that  are  at  constant  and  uniform  temperature  with  a  accuracy  of  perhaps  0.001  C  at 
temperatures  relevant  to  semiconductor  processing.  We  will  also  have  good  optical  access  to  the  samples. 
This  will  afford  the  opportunity  to  perfect  thermorefiectance  and  other  methods  of  non-invasive  temperature 
measurement.  Details  of  this  proposal  will  be  submitted  separately. 

3.  Effects  of  Temperature  and  Strain  on  Band  Offsets:  “The  Temperature  Dependence  of  Band  Offsets 
for  Semiconductor  Heterojunctions  in  General  and  for  the  Particular  Cases  of  AlAs-GaAs  and 
HgTe-CdTe,'  J.A.  Van  Vechten  and  K.J.  Malloy  (U.  of  Cal.,  Berkeley),  J.  Phys.:  Condens.  Matter. 
2,  281  (1990);  Thermal  Expansion  Contributions  to  Band  Gap  and  Band  Offset  Temperature 
Dependences."  K.J.  Malloy  (U.  of  NM)  and  J.A.  Van  Vechten,  J.  Vac.  Sci.  Technol.  (accepted)  TBP 
(1991). 

We  have  devised  a  simple  theory  for  the  temperature  and  strain  dependences  of  the  valence  and 
conduction  band  offsets  in  semiconductor  heterojunctions  using  the  thermodynamic  point  of  view.  We  find 
this  to  be  a  ma|or  advance  over  quantum  mechanical  approaches  to  the  same  problem.  We  have  offered 
a  resolution  to  conflicts  between  experimental  determinations  of  band  offsets  In  technically  Important 


6 


heterojunctions  by  reputable  groups  using  different  techniques  at  different  temperatures.  We  showed  that 
both  groups  can  be  right  for  the  conditions  they  applied.  We  have  also  developed  the  machinery  to  deal 
with  the  consequences  of  misfit  strain  and  dislocation  nucieation  on  the  electronic  properties  of 
heterojunctions  and  upon  the  point  defect  diffusion  near  them. 

4.  Elimination  of  DX  centers  from  AiGaAs  bv  Optically  Stimulated  Dislocation  Climb:  "Liquidation:  the 
Final  Solution  for  the  DX  Problem"  J.A.  Van  Vechten,  S.-B.  Kim,  and  Zhang  Q.-S.  and  R.  Matson 
(SERI)  submitted  to  the  16th  International  Conference  on  Defects  in  Semiconductors,  Lehigh,  PA 
7/91,  TBP. 

We  studied  the  effect  of  dislocation  climb  on  the  DX  center  concentration,  [DX],  in  n-type  AIGaAs. 
We  found  a  definite  decrease  in  [DX]  and  Increase  in  free  electron  concentration,  n,  concomitant  with  the 
climb,  which  was  optically  stimulated.  Clearly,  the  DX  centers  were  gettered  as  lattice  sites  were  added  to 
the  dislocation  structure  and  it  grew  dendritic  kinks  and  the  donor  dopant  atoms  (Te)  were  left  behind  with 
less  compensation.  We  hope  to  be  able  to  develop  this  into  a  practical  method  of  removing  the  DX  problem 
from  AiGaAs  and  other  wide  band  gap  compound  semiconductors. 

The  samples  used  in  this  study,  which  also  had  partial  support  from  the  NSF,  were  30%  Al  and 
grown  by  MOCVD  by  V.  Sundraman  of  Boeing  on  a  standard  <100>  GaAs  substrate.  Misfit  dislocations 
were  nucleated  by  scratching  the  surface  and  then  caused  to  glide  in  a  <011  >  direction  by  Kr+  laser 
scanning  irradiation  as  described  by  Monemar,  Van  Vechten  et  al.  (Phys.  Rev.  Lett.  11,  260  (1978)).  This 
optically  stimulated  glide  process  has  a  sharp  power  threshold  so  the  dislocations  were  stopped  at  boundary 
formed  by  an  Al  mask.  The  threading  parts  of  the  dislocation  networks  run  from  the  junction  up  3  microns 
to  the  sample  surface  near  the  mask  boundary.  We  caused  these  threading  dislocations  to  climb  by 
illuminating  the  sample  with  a  laser  or  with  white  light  at  intensities  below  the  glide  threshold,  as  also 
described  previously.  We  determined  [DX]  both  by  DLTS  and  by  capacitance  transient  spectroscopy  at 
constant  temperature,  C(t),  and  determined  n  by  C-V  profiling.  In  order  to  have  adequate  junction  area  for 
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these  electrical  measurements,  we  made  the  mask  holes  0.2  mm  x  0.2mm,  which  is  far  larger  than  the  2 
micron  size  envisioned  for  practical  gettering  of  discrete  devices.  Consequently,  the  gettering  process  took 
about  10  hours,  a  factor  of  10,000  more  than  for  the  practical  masks.  The  degree  of  getting  was  only  about 
50%  for  these  large  mask  samples  but  it  is  clear  that  our  initial  attempts  are  far  from  optimum.  We  hope 
to  have  much  more  complete  gettering  with  the  smaller  masks  and  with  improved  processing. 

Lang  et  al.  (Phys.  Rev.  Lett.  42  1353  (1979)  had  previously  reported  that  the  climb  of  a  threading 
dislocation  in  AIGaAs  reduced  [DX],  but  this  report  seems  to  have  been  generally  forgotten  in  recent  years. 
Ours  seems  to  be  the  first  attempt  to  combine  these  two  reports  from  1 978-1 979  to  affect  defect  engineering 
by  optically  induced  dislocation  glide  and  climb.  Part  of  the  reason  that  Lang  et  al.’s  report  has  fallen  from 
notice  is  that  many  workers  in  the  field  have  believed  models  for  the  DX  center,  such  as  that  of  Chadi  and 
Chang,  which  hold  that  the  deep  center  and  the  donor  atom  are  one  and  the  same.  If  this  *No  X'  hypothesis 
were  true,  it  would  not  be  possible  to  getter  the  deep  center  without  reducing  the  net  donor  doping,  n,  also. 
We  submit  that  this  present  result,  along  with  several  others,  completely  disprove  all  “No  X“  theories  for  the 
DX  center  in  AIGaAs. 

5-  Proper  Analysis  of  the  Capacitance  Transient  due  to  Multiple  Carrier  Emission  from  a  Highly 
Compensating  Deep  Level,  e.o..  DX:  “Capacitance  Transient  Due  to  Multiple  Majority  Carrier 
Emission  from  DX  Centers  in  AIGaAs,"  S.-B.  Kim,  J.A.  Van  Vechten,  and  Zhang  Q.-S.,  TBP. 

Because  of  the  controversial  nature  of  our  result  for  the  effect  of  dislocation  climb  upon  DX  center 
concentration,  [DX],  just  described  above,  we  made  a  careful  study  of  the  capacitance  transient  at  constant 
temperature,  C(t),  technique  which  we  used  to  determine  [DX].  Almost  everyone  who  uses  this  technique 
applies  an  analysis  which  assumes  (see  e.g.,  D.  K.  Schroder  "Semiconductor  material  and  device 
characterization*  J.  Wiley,  New  York,  1990,  Chapt.  7,  especially  pp.  306)  that  the  deep  levels  trap  and  emit 
free  carriers  in  one  transition  only,  that  the  applied  voltage  in  the  reversed  biased  condition,  VA,  is  much 
larger  than  the  built  in  voltage,  Vbi,  of  the  junction  employed,  and  that  the  density  of  deep  centers  is  much 
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less  than  the  net  doping  level.  It  is  now  firmly  established  that  the  DX  center  traps  at  least  2  and  more  likely 
3  electrons  (Van  Vechten,  J.  Phys.  Cond.  Matter.  1,  5171  (1989);  Gibart  et  al.,  Phys.  Rev.  Lett.  65,  1144 
(1990)).  Our  C(t)  experiments  and  those  of  many  other  groups  are  done  with  -  VA  ~  Vbj.  The  DX  centers 
compensate  typically  2/3  or  more  of  the  donor  doping,  which  is  why  there  is  so  much  interest  in  the  DX 
center.  Thus,  all  three  of  these  assumptions  made  in  the  standard  analysis  are  Invalid  for  the  typical 
experiment  to  determine  [DX]  in  AIGaAs. 

We  revised  the  analysis  of  C(t)  to  take  account  of  the  the  facts  that  the  compensation  ratio  can  be 
large,  *  VA  can  be  about  the  same  as  Vbi,  and  that  multiple  carriers  can  be  trapped  and  emitted  in  several 
different  transitions.  We  treated  the  depletion  layer,  W,  In  two  sections.  The  first,  WFP,  Is  that  which  is 
depleted  by  Vbi  alone  when  VA  =  0,  as  during  the  Tilling  pulse’.  In  WFP  we  assume  that  the  compensation 
ratio,  i  e.,  1  minus  the  ratio  of  the  Ionized  dopant  concentration  to  the  net  charge  density,  is  Invariant  in  time 
and,  for  a  homogeneous  sample,  constant  in  space.  We  call  this  compensation  ratio  g.  When  the  junction 
is  reverse  biased,  the  depletion  width  is  wider  and  the  region  beyond  WFP  contains  the  majority  carriers  that 
have  been  trapped  by  the  deep  centers.  In  our  case  these  are  the  electrons  that  have  been  trapped  by  the 
DX  centers.  Thus,  the  compensation  ratio  in  this  expanded  part  of  the  depletion  layer  is  generally  larger  than 
that  in  Wpp  and  varies  in  time  as  the  trapped  carriers  (electrons)  are  emitted.  We  call  this  compensation 
ratio  f(t)  and  take  it  to  be  uniform  across  the  expanded  part  at  any  given  time.  With  increasing  time  under 
reverse  bias,  f(t)  approaches  g.  As  this  happens,  W,  which  we  determine  by  double  integrating  the  Poisson 
integral  of  these  two  charge  layers,  decreases  and  the  junction  capacitance,  C,(t),  increases.  It  is  Cj(t)  that 
is  actually  measured. 

We  have  derived  generally  expressions  for  Cj(t)  for  this  two  layer  model  of  the  depletion  region  for 
the  cases  that  a)  2  electrons  are  emitted  in  one  ’negative-U"  transition  (as  per  the  Chadi  Chang  model);  b) 
2  electrons  are  emitted  in  two  correlated  transitions;  c)  3  electrons  are  emitted  three  correlated  transitions 
(as  per  the  Lang-Logan-Jaros,  Van  Vechten  model)  d)  3  electrons  are  emitted  in  two  transitions,  one  of 
which  is  negatlve-U;  and  e)  3  electrons  are  emitted  In  three  independent  transitions.  These  expressions  are 
somewhat  complicated  by  the  standards  of  the  era  before  the  widespread  availability  of  micro  computers, 
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which  is  probably  why  this  analysis  was  not  done  previously.  Howevei,  we  have  developed  and  are 
distributing  a  robust  fitting  program  that  analyses  raw  data  in  terms  of  any  of  these  assumptions  one 
chooses  to  use.  This  progiam  runs  on  a  standard  PC  (with  math  coprocessor)  in  a  few  seconds. 

When  applied  to  cur  own  C(t)  data  for  the  samples  described  above,  we  find  that  for  one  fixed 
temperature  (e.g.,  131  K)  all  five  of  the  above  assumptions  can  be  parameterized  to  fit  the  data  well  enough 
to  pass  a  95%  statistical  confidence  test.  However,  by  Arrhenius  analysis  of  the  computer  fitted  parameters 
through  the  temperature  range  125  to  137  K,  we  can  make  a  strong  argument  that  options  d  and  e  should 
be  discarded  for  the  DX  center.  Furthermore,  the  conclusion  that  [DX]  has  been  reduced  by  the  dislocation 
climb,  and  how  much  this  reductior  is  for  a  given  sample,  is  quite  Insensitive  to  which  of  the  frve  choices 
one  makes.  Thus,  we  are  confident  r  the  validity  of  the  conclusion  of  the  previous  section  that  optically 
stimulated  dislocation  climb  does  reduce  [DX]  while  simultaneously  increasing  n. 

6.  C-V  Analysis  and  SPICE  Modeling  of  ACTFEL  Devices:  "Capacitance-Voltage  Characteristics  of 
Alternating-Current  Thin-Film  Electron  Luminescent  Devices,"  R.C.  McArthur,  J.D.  Davidson,  J.F. 
Wager,  I.  Khormaei,  and  C.N.  King,  Appl.  Phys.  Lett.  §6.  (1990),  "Characterization  of  ZnS:Mn 

AC  Thin-Film  Electroluminescent  Devices  by  Capacitance-Voltage  Analysis,"  R.C.  McArthur,  J.D. 
Davidson,  I.  Khormaei,  J.F.  Wager,  and  C.N.  King,  Acta  Poly  Technica  Scandinavia  Ph  170.  181 
(1990);  "Aging  Instabilities  of  ZnS:Mn  AC  Thin-Film  Electroluminescent  Devices,"  J.D.  Davidson,  J.F. 
Wager,  I.  Khormaei,  and  C.N.  King,  Acta  Polytechnica  Scandinavia  Ph  170.  185  (1990);  "Electiical 
Characterization  and  Modeling  of  Alternating-Current  Thin-Film  Electroluminescent  Devices,"  J.D. 
Davidson,  J.F.  Wager,  I.  Khormaei,  C.N.  King,  and  D.  Williams,  IEEE  Trans.  Electron  Devices 
(submitted) 

The  standard  method  for  characterizing  the  electrical  properties  of  alternating-current  thin-f.’m 
electroluminescent  (ACTFEL)  display  devices  is  *v  charge-voltage  (Q-V)  analysis.  We  have  developed  a 
complementary  technique,  the  capacitance-voliage  (C-V)  technique  for  ACTFEL  characterization.  Since 
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ACTFEL  device  operation  depends  upon  the  application  of  a  large  amplitude  (300-450  V  peak-to-peak)  AC 
waveform,  conventional  small-signal  C-V  analysis  is  precluded.  Instead,  we  perform  our  C-V  analysis  by 
sampling  the  instantaneous  voltage  across  a  current -limiting  resistor  placed  in  series  with  the  ACTFEL  device 
and  processing  these  sampled  signals  digitally  to  obtain  the  C-V  curve.  Q-V,  Q-t,  and  l-t  analysis  can  be 
performed  in  the  same  manner. 

Our  motivation  for  developing  this  technique  is  to  study  the  properties  of  point  defects  in  the  ll-VI 
semiconductor  ZnS,  as  manifest  by  the  electrical  properties  of  commercial  devices.  Specially  fabricated, 
commercial  quality  devices  are  fabricated  by  Planar  Systems,  Inc.,  Beaverton,  Oregon,  for  this  purpose. 

The  slope  of  the  C-V  curve  in  the  transition  region  is  shown  to  be  related  to  the  density  of  interface 
states  at  the  ZnS-insula*or  interface.  Commercial  devices  manufactured  by  Planar  Systems  are  found  to 
have  moderate  interface  state  densities  which  are  independent  of  the  operatic  time  of  the  display.  In 
contrast,  commercial  devices  supplied  by  Finlux,  Helsinki,  Finland,  initially  have  a  low  density  but  the 
interface  state  density  increases  dramatically  as  the  devices  are  operated.  We  are  currently  trying  to 
understand  the  physical  nature  of  these  interface  defects  and  their  dynamics  with  operation  time. 

We  have  also  recently  applied  SPICE  modeling  in  conjunction  with  a  simple  equivalent  circuit  for 
the  ACTFEL  device  in  order  to  assist  in  the  interpretation  of  C-V  curves.  We  have  found  SPICE  modeling 
to  be  invaluable  as  a  tool  for  clarifying  ACTFEL  device  physics  issues. 

The  C-V  technique  has  proven  to  be  useful  for  experimentally  investigating  Lie  aging  instabilities 
inherent  in  present  art  ACTFEL  devices.  We  are  performing  variable-temperature  aging  studies  on  ZnS:Mn 
ACTFEL  devices  in  an  attempt  to  deduce  the  microscopic  mechanism  responsible  for  aging  instabilities. 
Preliminary  analysis  of  the  kinetic  data  suggests  that  two  aging  mechanisms  are  operational.  One 
mechanism,  which  dominates  at  low-temperatures  and  has  a  small  activation  energy,  we  attribute  to  electron 
injection  and  trapping  a»  the  ZnS/insulator  interface.  The  seconc  mechanism  is  important  at  and  above 
room  temperature  and  is  characterized  by  a  larger  activation  energy;  we  attribute  this  aging  mechanism  to 
atomic  migration  or  intei diffusion  which  is  driven  by  hot  electron  recombination. 
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We  are  applying  atomistic  thermodynamic  modeling  to  these  problems  in  order  to  attempt  to 
understand  on  a  microscopic  scale  the  nature  of  these  device  properties. 

7.  DX  in  AIGaAs:  Transient  Decay  of  Persistent  Photoconductivity  in  AI^GaQ  jAs‘  T.W.  Dobson, 

LV.A.  Scalvi,  and  J.F.  Wager,  J.  Appl.  Phys.  S2.  601(1990). 

One  of  the  peculiarities  of  the  DX  center  In  AIGaAs  is  the  fact  that  it  exhibits  persistent 
photoconductivity  (PPC)  at  low  temperatures  (s  50K).  At  Intermediate  temperatures  (80-1 00K)  the  PPC  is 
not  entirely  persistent  and  exhibits  a  decay  with  time.  We  have  measured  this  transient  decay  of  PPC 
(TDPPC)  as  a  means  of  exploring  the  donor-  or  acceptor-like  nature  of  the  DX  center.  We  analyze  our 
experimental  data  by  performing  a  computer  simulation  which  Involves  the  nonequPibrium  capture  kinetics 
of  electrons  from  the  conduction  band  into  the  DX  center  and  a  concomitant  change  in  mobility  due  to 
ionized  impurity  scattering.  We  then  simulate  six  models  for  the  DX  center.  We  find  excellent  quantitative 
agreement  between  experiment  and  the  simulation  for  the  DX  model  of  Chadl  and  Chang  [Phys.  Rev.  B29, 
1 0063  (1989)],  which  assumes  DX  to  be  a  negative-U  acceptor  with  intrinsic  self-compensation,  plus  a  photo- 
induced  shallow  donor. 

8.  DX  in  AIGaAs:  S.B.  Kim  and  J.F.  Wager,  J.  Appl.  Phys.  (in  preparation). 

The  low-temperature  minority  carrier  capture  cross  section  for  the  DX  center  has  been  measured 
for  the  first  time.  The  minority  carrier  capture  (MCC)  experiment  is  accomplished  by  photo-injection  of  holes 
into  the  space  charge  region  of  a  semitransparent  Schottky  barrier  on  AIGaAs.  We  find  the  capture  cross 
section  to  be  of  magnitude  2x1  O'1 9  cm2  and  Independent  of  temperature  over  the  range  30-70K.  According 
to  the  theory  of  Ridley  [Solid  State  Electronics  21.  1319  (1978),  J.  Phys.  C.  14,  1255  (1981)]  and  the 
experimental  work  of  Henry  and  Lang  [Phys.  Rev.  B15, 989  (1977)],  a  low-temperature  capture  cross  section 
of  this  magnitude  and  which  is  temperature-independent  is  Indicative  of  hole  capture  by  a  coulombically 
neutral  trap.  Thus,  we  conclude  that  hole  capture  occurs  into  a  neutral  DX  trap  state.  We  note  that  this 
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result  can  be  reconciled  with  the  DX  model  of  Chadi  and  Chang  [Phys.  Rev.  B22,  10063  (1989)]  If  the 
capture  is  postulated  to  occur  by  the  f  Jlowing  two-step  process  in  which  the  second  step  is  rate-limiting: 
Stepi:  h+  +  DX'  =  DX° 

Step  2:  h+  +  DX°  =  d+  (rate-limiting) 

9.  Statistical  Thermodynamics  of  Ballistic  Hooping:  "A  Statistical  Thermodynamic  Derivation  of  the 
Ballistic  Model  for  Vacancy  Migration,"  J.F.  Wager,  Phil.  Mag.  A  (1991),  in  press. 

Much  of  our  work  on  defect  dynamics  is  based  on  the  validity  of  the  Ballistic  Model  (BM)  [J.A.  Van 
Vechten,  Phys.  Rev.  12,  1247  (1975)]  which  asserts  (essentially)  that  the  enthalpy  for  vacancy  migration  Is 
determined  by  the  kinetic  energy  of  the  atom  which  hops  into  the  vacancy.  We  have  previously  extended 
the  BM  hypothesis  to  calculate  the  entropy  of  migration  for  nearest-neighbor  hopping  [T.W.  Dobson,  J.F. 
Wager,  and  J.A.  Van  Vechten,  Phys.  Rev.  B12.  2962  (1989)]. 

The  primary  impetus  for  this  work  is  that  the  migration  entropy  is  estimated  by  our  previous  method 
[Dobson  et  al.]  is  anomalously  large  when  compared  to  known  trends  for  self-diffusion  in  elemental  metals. 
The  statistical  thermodynamic  derivation  leads  to  a  migration  enthalpy  and  entropy  qualitatively  similar,  but 
not  Identical,  to  that  given  previously.  The  statistical  thermodynamic  correction  to  the  migration  enthalpy 
i«  of  negligible  importance.  In  contrast,  the  correction  to  the  entropy  of  migration  which  arises  from  the 
statistical  thermodynamic  derivation  is  found  to  be  significant  and  must  be  employed  in  order  to  obtain 
agreement  with  experimental  values  for  eiementai  metals. 

10.  Energetics  of  Self-Diffusion  in  GaAs:  "Energetics  of  Self-Diffusion  in  GaAs,"  J.F.  Wager,  J.  Appl. 
Phys.  §2,  3022  (1991). 

In  order  to  provide  insight  into  the  nature  of  the  mechanism  of  self -diffusion  in  GaAs,  two  kinds  of 
atomistic  thermodynamic  calculations  are  performed.  First,  the  activation  enthalpy  of  the  saddle  point 
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configuration  for  various  modes  of  vacancy-mediated  self-diffusion  is  estimated.  The  primary  conclusion 
of  the  analysis  is  that  if  vacancies  are  available  in  equilibrium  concentrations  vacancy  migration  by  second 
nearest-neighbor  hopping  (2nnh)  is  found  to  be  the  energetically  most  favorable  mechanism  of  self-diffusion, 
independent  of  the  Fermi-level  position.  The  second  type  of  atomistic  thermodynamic  calculation  performed 
is  estimation  of  the  activation  entropy  for  various  modes  of  vacancy  migration.  The  primary  conclusion  of 
the  analysis  is  that  vacancy  migration  by  2nnh  is  consistent  with  activation  entropies  of  order  6  kB  whereas 
vacancy  migration  by  nearest-neighbor  hopping  (nnh)  is  characterized  by  much  larger  activation  entropies 
of  order  20  kB. 

Our  activation  enthalpy  and  entropy  analysis  leads  to  the  conclusion  that  experimental  self-diffusion 
and  Interdiffusion  data  with  abnormally  large  diffusivity  prefactors  and  activation  energies  [as  summarized, 
for  example,  by  T.Y.  Tan  and  U.  Gdsele,  Appl.  Phys.  Lett  52,  1240  (1988)]  cannot  be  explained  in  terms  of 
any  vacancy  mediated  self-diffusion  mechanism  in  which  vacancies  are  present  in  equilibrium  concentra¬ 
tions.  Rather,  we  propose  that  Ga  self-diffusion  or  interdiffusion  processes  which  are  characterized  by  large 
prefactors  (~107-108  cm2  s'1)  and  activation  energies  (~6  eV)  involve  cation  Frenkel  pair  formation.  This 
mode  of  self-diffusion  is  expected  to  be  dominant  when  surface  vacancy  generation  is  precluded  and 
vacancies  can  only  be  generated  In  the  bulk  of  the  semiconductor  by  the  creation  of  an  interstitial-vacancy 
pair. 

The  results  of  our  activity  are  discussed  in  detail  in  the  attached  reprints  and  preprints. 
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